Ray model and ray-wave correspondence in coupled optical microdisks 
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We introduce a ray model for coupled optical microdisks, in which we select coupling-efficient 
rays among the splitting rays. We investigate the resulting phase-space structure and report island 
structures arising from the ray-coupling between the two microdisks. We find the microdisks's re- 
fractive index to influence the phase-space structure and calculate the stability and decay rates of 
the islands. Turning to ray-wave correspondence, we find many resonances to be directly related 
to the presence of these islands. We study the relation between the (ray-picture originating) is- 
land structures and the (wave-picture originating) spectral properties of resonances, especially the 
leakiness of the resonances which is represented as the imaginary part of the complex wave vector. 

PACS numbers: 42.55.Sa, 05.45.Mt, 42.25.-p, 42.60.Da 



I. INTRODUCTION 

With the advances in material science and nano- 
fabrication techniques, coupled optical microdisks have 
recently attracted much interest especially in the context 
of device applications such as photonic molecules and 
coupled-resonator optical waveguide (CROW). The ter- 
minology 'photonic molecule' is introduced as an optical 
analog to chemical molecules. Their molecule-like struc- 
ture arises from pairs of interacting optical microdisks 
(that, in the analogy, represent the individual atoms) 
that are coupled by narrow channels [3-Q . CROW is the 
chain version of an optical molecule, and this chain of 
coupled optical microdisks can be used as an optical de- 
lay line through weak coupling between localized optical 
high-Q cavities 0, HJ. In addition, coupled optical mi- 
crodisks have also been studied for theoretical prospects 
and experimental realizations in the contexts of optical 
mode coupling and producing directional light emission 



The ray picture and the associated ray dynamics have 
been investigated extensively in optical microcavities 
since asymmetric resonant cavities, whispering-gallery 
resonators with smooth deformations of a circular bound- 
ary shape, were introduced in order to break the rota- 
tional invariance of disk resonators and obtain directed 
emissions from optical microcavities [3, [Hi. In those 
slightly deformed microcavities, light is predominantly 
emitted from the boundary points of highest curvature 
and the direction is tangential to the boundary. From a 
ray dynamical viewpoint, this is a result of tunneling of 
the rays supporting the whispering gallery mode (WGM) 
confined in the Kolmogorov-Arnold-Moser (KAM) tori 
through the lowest dynamical barrier. As the microcav- 
ity is strongly deformed, the ray dynamics becomes more 
chaotic and the direction of light emission is determined 
by the unstable manifold structure fl6l - [2ll ]. There are 
many recent successes in understanding and designing 
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microcavities based on a ray dynamical analysis that ex- 
plained the existence of high-Q modes with directional 
emission characteristics. 

However, in coupled optical microcavities, the identi- 
fication of a ray dynamical model and consequently the 
establishment of ray-wave correspondence are not simple 
because of the intrinsic refractive ray splitting. Quantum 
chaos in systems with ray splitting has been studied by 
ray models with stochastic selection rule [2^ - |28| and for 
special resonator geometries such as annular cavities [2!| . 
In this paper, we introduce the ray model for coupled op- 
tical microdisks, in which we select coupling-efficient rays 
among the splitting rays by imposing a deterministic se- 
lection rule and study ray-wave correspondence based on 
this ray model. In Sec. II we explicate our ray model in 
coupled dielectric disks and show how island structures 
in phase space originate from the coupling between two 
disks. We also obtain the stabilities and decay rates of 
principal periodic orbits as well as the island structures. 
Regular modes localized on islands are presented in Sec. 
Ill and we study the relation between classical structures 
and resonance modes systematically. Finally, we summa- 
rize the results in Sec. IV. 



II. RAY MODEL WITH DETERMINISTIC 
SELECTION RULE FOR COUPLED 
DIELECTRIC DISKS 

The ray model of a dielectric disk is very simple. If the 
incident angle of a ray is larger than the critical angle for 
total internal reflection, the ray totally reflects and, due 
to the rotational symmetry and the related conservation 
of angular momentum, circulates inside the disk forever 
with the same incident angle. If the incident angle of a 
ray is smaller than the critical angle, the ray reflects and 
transmits partly - the ray is refractively splitted. The 
loss of rays to the outside implies that the overall inten- 
sity inside the disk will be reduced according to (a gen- 
eralized) Fresnel's law [3(|, but apart from this, all rays 
inside the disk will circulate keeping the initial angle as 
the incident angle for each reflection and no additional 
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FIG. 1: (color online) (a) Coupled dielectric disks. Simple pe- 
riodic orbits (red lines) in coupled dielectric disks, (b) Verti- 
cal bouncing ball type periodic orbit, (c) Horizontal bouncing 
ball type periodic orbit, (d) Hexagonal-shaped periodic orbit. 
The blue dotted lines represent the associated quasiperiodic 
(or periodic) trajectories in a single disk. 



features arise. The ray model of a two-disk billiard con- 
sisting of two closed disks is also very easy to understand: 
If a ray is inside one of the disks, it will again circulate 
inside this disk with the same incident angle. If a ray 
approaches the disks from outside, it will escape after 
one scattering at the (hard-wall) disks. Note that there 
exists one unstable periodic orbit in the system, namely 
the shortest path between the two disks. 

However, the ray model of two coupled dielectric (i.e., 
partially open) disks with radii Rl and Rr, refractive 
indices til and Hr, and a separation d [cf. Fig. [1] (a)] is 
not simple because the ray which emits from one disk 
can, in principle, enter the other disk. At the recurrence 
of the ray into the first disk, the incident angle will be 
varied, and the dynamical properties of a ray model for 
coupled dielectric disks become plentiful, as will become 
visible in a more complicated phase-space structure. An 
additional complication arises from the fact that the ray 
splitting leads to an exponential multiplication of rays, 
namely to 2™ trajectories after n bounces, that are, in 
practice, impossible to follow individually. Therefore, a 
selection of rays has to be made, based on a selection rule 
that picks those rays that determine the system dynamics 
and properties. There are many possibilities for imposing 
a selection rule; we have chosen a deterministic selection 
rule that emphasizes the coupling between the two disks: 
Among the two splitted rays, we select the one that is 
more effective for the coupling between the two disks - 
if one of the splitted rays (e.g. the transmitted part of 
a ray approaching the disk boundary from inside one of 
the disks) will reach the other disk, we choose this one. 
If this ray would, however, escape the system, we choose 
the ray that remains inside (enters into) one of the disks. 
Totally reflected rays do not exhibit ray splitting and 
remain in the respective disk forever. 

A more precise formulation of our deterministic selec- 
tion rule reads as follows: If \p\, modulus of the sine of 
the incident angle, of an initial ray inside the left disk 
is larger than the critical p c = no /til = l/ n f° r total 
internal reflection, the ray circulates inside the left disk 



and does not reach the right disk. If the initial ray is lo- 
cated in the open (refractive) region, i.e. — p c < p < p c , 
the ray reflects and transmits partly according to Fres- 
nel's equation [3(1 HH ■ The reflected ray remains inside 
the left disk but the transmitted ray either escapes from 
the system or approaches the right disk. We select the 
reflected ray if the transmitted ray emits from the sys- 
tem but the transmitted ray if it reaches the right disk. 
The ray that hits the right disk also reflects and trans- 
mits partly and here we always discard the reflected ray 
and keep the transmitted ray that always enters into the 
second disk. We repeat this selective process whenever 
the rays meet the dielectric boundary. Note that, con- 
sequently, none of the kept light rays ever escapes the 
system. Contrary to the previous studies on ray splitting 
model with a stochastic selection rule [22l - [28j . we intro- 
duce here a ray model with a deterministic selection rule 
(RMDS) and apply it throughout this paper. 

We first consider periodic orbits in coupled dielectric 
disks in the framework of the RMDS ray model. In com- 
parison to the single disk, coupling between two disks due 
to openness produces many new periodic orbits. FigurcQ] 
(b)-(d) show three types of simple periodic orbits in cou- 
pled dielectric disks. The vertical bouncing ball type pe- 
riodic orbits of Fig. [1] (b) are the same as those of a single 
dielectric disk, except for a rotational symmetry break- 
ing resulting from the coupling. The horizontal bounc- 
ing ball type periodic orbit of Fig. [T] (c) consists of the 
bouncing ball type periodic orbit in the individual disks 
and the shortest path between them which is the only 
periodic orbit outside the two disks. As these periodic 
orbits [red lines of Fig. [T] (b) and (c)] are made by com- 
bining periodic orbits which exist in disks with a closed 
boundary condition, they always exist in the coupled di- 
electric disk system, independent of system parameters. 
The hexagonal-shaped periodic orbit of Fig.[T](d) is a new 
periodic orbit for which no corresponding periodic orbits 
exist in a closed boundary condition. Consequently, its 
existence and (geometric) properties will depend on the 
system parameters. 

We now study the ray dynamics in our RMDS-model 
in terms of the Poincare surface of section (PSOS) that 
is obtained by plotting the position s (arclength along 
disk boundary, cf. Fig. Q] (a)) and the correspond- 
ing p, sine of the angle of incidence, for each reflection 
point at the boundary of the left disk. The PSOS for the 
symmetric system (til/iio = ur/uq = n and Rl = Rr, 
all lengths are measured in units of the interdisk spac- 
ing d = 1.0) is shown in Fig. [5J It is obtained from 
initial rays uniformly distributed in phase space and the 
ray trajectories resulting from the RMDS for coupled di- 
electric disks. In the region of total internal reflection 
(\p\ > p c ), the ray dynamics follows horizontal lines that 
are single-disk WGMs. Differences from the single-disk 
behavior occur, as expected, in the open region (\p\ < p c ) 
of phase space where the dielectric coupling between the 
two disks becomes effective. We emphasize that the ray 
dynamics will go through a transient behavior before the 
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FIG. 2: (color online) (a) The ray attractors in the left disk 
of coupled dielectric disks with interdisk distance d = 1.0, 
ratio of radii R = Rl/Rr = 1.0, and refractive indices 
n — n^/no = ur/uo = 2.0. The blue dashed line marks 
the critical line p c = 1/n = 0.5 for total internal reflection. 
Note that the right and left part of the attractor are related 
by the ray dynamics of the ORM within the left disk, (b) En- 
larged ray attractor corresponding to the red square of (a). 
The arrows A (lower red), B (upper green), and C (middle 
blue) mark the attractors that correspond to the red, green, 
and blue trajectories, respectively, shown on the right. 



structure displayed in Fig. [2] is reached in the stationary 
regime. Since the ray dynamics within the RMDS-model 
depends on the system parameters (geometrical details, 
refractive index ratios), so will the details of the PSOS 
in this region, see also Fig. [3] and the discussion below. 
For the symmetric geometry, the PSOS structure in the 
open region forms an attractor in the stationary regime, 
see Fig. [2] (a), i.e., the rays with initial \p\ smaller than 
p c go to the attractor after a transient time. The black 
areas in Fig. [2] (a) represent the attractor which actually 
has regular island structures as shown in Fig. [2] (b). For 
instance, rays started at positions (s,p) = (0.0,0.15) and 
(0.0,0.36) on the attractor form the red (A) and green 
(B) trajectories in configuration space (upper and cen- 
tral panel on the right of Fig. [2]) and the corresponding 
sets in phase space (marked by arrows). A ray with ini- 
tial position (s,p) = (0.0,0.45) enters the blue (C) set 
(lower panel on the right of Fig. [2]) after the ray travels 
the open region for the transient time and then stays on 
the attractor forever. 

Figure [3] shows the attractor structure in the PSOS 
for different refractive indices and the same symmetric 
geometry. Clearly, the details of the attractor depend 
crucially on the underlying ray dynamics, that is, (be- 
sides a geometry dependence) the index n of refraction. 
For n = 1.5, cf. Fig. 02a), the attractor covers a broad 
range of phase space and its fractal structure represents 
the fully chaotic ray dynamics. The onset of island for- 
mation can be seen for n = 1.8 in panel (b). Further 
increasing n enlarges the islands and adds to the richness 
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FIG. 3: PSOSs in the open region (\p\ < p c ) on the phase 
space when (a) n — 1.5, (b) n = 1.8, (c) n = 2.2, and (d) 
n = 3.0. 



of their structure, recall the detail shown in Fig. &b) for 
n = 2. This behavior is reminiscent of the evolution of 
KAM islands in the context of chaotic dynamics, e.g., 
as an external parameter such as the kicking strength in 
the kicked rotator is changed. In the present case, it is in 
particular the n-dependence of the selection rule for the 
RMDS that effectively changes the underlying (RMDS) 
map. The detailed description of the attractors will be 
reported elsewhere. 

The attracting behavior of the RMDS is caused by the 
underlying deterministic selection rule, as illustrated in 
Fig- El The attractor structure depends on the parameter 
n that in turn determines the selection rule. The basin 
of attraction is formed by the open region in phase space 
that extends between momenta p = ±l/n (rays started 
outside this region will remain above/below the critical 
line and inside the disk they started in forever). When 
a ray reaches the attractor after some transition time, it 
then has to stay on the attractor forever, according to the 
rules of the RMDS: Rays that would leave the attractor 
are never selected. Note that this also implies that no 
attractors are formed if the full ray dynamics, combining 
the features of the RMDS and the ordinary ray model 
(ORM) of the individual disks, is taken into account. 
Therefore, the attractors indicate invariant sets to which 
ray trajectories approach in the dynamics determined by 
the RMDS. Notice that these attractors are similar to the 
"quasi-attractors" in a piecewise smooth area-preserving 
map for which noninvertibility was found to induce the 
phase space collapse [H, [H| . We also point out that the 
attractors that we study here are different from those 
described in Ref. (HJ which resulted in the context of 
optical microcavities with a corrected ray dynamics. 

A striking feature of the RMDS attractors is the exis- 
tence of regular islands. For instance, although the pe- 
riodic orbit of Fig. [1] (c) consists of marginally stable 
periodic orbits inside the disks and an unstable periodic 
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orbit connecting the disks, the combined RMDS periodic 
orbit is marginally stable with elliptic structure as shown 
in Fig. [2] (b). This is confirmed when studying its sta- 
bility via its monodromy matrix M that, in terms of the 
Birkhoff coordinates, is defined by [36 1 



dsj 
5pj 



M 



5sj 
S Pj 



ran mi2 
m 2 \ m 2 2 



(1) 

where Sj is the arclength s at the jth bounce point and 
angular momentum pj is p at the jth bounce point. A 
trajectory starts from a position Sj-% with momentum 
Pj-i and ends at a position Sj with momentum pj. We 
then take another trajectory starting with slightly differ- 
ent initial conditions Sj-i +6sj-\ andpj-i +Spj-\. The 
deviations at the end position are obtained from the de- 
viations at the starting position in linear approximation 

as E q . mj mm- 

To obtain the monodromy matrix of the periodic or- 
bit of Fig.[T](c), we divide this orbit into three elements, 
namely the trajectories inside and in between the disks, 
and the refraction at the dielectric interfaces [39M4l| . 
First, we obtain the monodromy matrix of the trajec- 
tory of the periodic orbit inside a disk. The monodromy 
matrix for a left (right) disk is thus 



M 



L(R) 



2R 



L(R) 



(2) 



From the relation Tr M L ( R ) = 2, we confirm that the 
trajectory inside the disks are marginally stable. Sec- 
ondly, we obtain the monodromy matrix of the trajec- 
tory in between the two disks. After some algebra, the 
corresponding monodromy matrix reads 



Mi = - 



1 



i 

Rl 



,1 

Rl 



RlRl 



l+if 

Rr 



(3) 



Since Tr Mi < — 2 if d > 0, this part of the trajectory is 
always unstable. Note that up to this point the problem 
is equivalent to the corresponding orbit in a two-disk bil- 
liard 42j . In order to obtain the full monodromy matrix 
for our couplcd-dielectric-disks problem, we need in addi- 
tion the monodromy matrices at the dielectric interfaces 
where the light is refracted. Since, according to Sncll's 
law, 6pj = nSpj^i if a ray goes from inside to outside a 
disk and Spj = ^Spj-i in the opposite case, the respec- 
tive monodromy matrices are 



and 



Mi 



M t 



1 

n 



1 
i 



(4) 



(5) 



Finally, the monodromy matrix of the horizontal bounc- 
ing ball type periodic orbit of Fig. [1] (c) is 

M = M l Mb 2 MiMb 1 M r M r Mb 2 MiM Bi M l . (6) 
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FIG. 4: (color online) (a) Contour plots which satisfy the con- 
ditions of Tr M = 2 (black line) and Tr M = -2 (red dashed 
lines) on (Rr, n)-parameter space when d = 1.0. PSOSs at 
(b) n = 1.98 (region I), (c) n = 2.5 (region II), and (d) n = 3.0 
(region III) when d — 1.0 and Rr = 1.2. 



Although Mbj and Mb 2 induce area-expanding and 
area-contracting properties, respectively, upon transmis- 
sion of rays through dielectric interfaces, the final lin- 
ear mapping for the full periodic orbit satisfies the area- 
preserving property with detM = 1. 

The stability of the horizontal periodic orbit depends 
on the value of Tr M. In the case of |Tr M| > 2, 
the periodic orbit is linearly unstable but in the case of 
|Tr M| < 2, the periodic orbit is linearly stable. In order 
to classify the linear stability of this orbit, contour plots 
satisfying the condition Tr M = ±2 in (Rr, n)-parameter 
space (with d = 1.0) are shown in Fig. 0] (a). The lin- 
ear stability can be classified into three different types 

mull- 
in region I with Tr M > 2, the two eigenvalues of M 
are a pair of reciprocal positive real values and the regular 
(ordinary) hyperbolic periodic orbit is linearly unstable. 
In region II with — 2 < Tr M < 2, the two eigenvalues of 
M form a complex conjugate pair on the unit circle and 
the corresponding elliptic periodic orbit is linearly stable. 
In region III with Tr M < —2, the two eigenvalues of 
M are a pair of reciprocal negative real values and the 
inverse (reflection) hyperbolic periodic orbit is linearly 
unstable. The existence of the periodic orbit, (s,p) = 
(0,0), is independent of the system parameters but the 
stability of the periodic orbit is determined by the system 
parameters in coupled dielectric disks. 

To confirm the stability diagram of Fig. [4] (a), we 
obtained numerically the PSOSs in the close vicinity 
of the periodic orbit, more precisely around its point 
(s,p) = (0,0) at the left disk boundary. Figure 0] (b), 
(c), and (d) show the PSOSs in the different stability re- 
gions of Fig. 2] (a). We start in region I with a suitable 
parameter n — 1.98 such that the horizontal periodic or- 
bit possesses a hyperbolic structure. As n decreases, the 
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scparatrix is broken and then the island structure near 
(0,±0.15) become smaller and finally disappear. In re- 
gion II, the periodic orbit has elliptic character. In region 
III, the periodic orbit has again a hyperbolic structure. 
As Rr increases, the separatrix is also broken and then 
the island structure near (±0.05,0) become smaller and 
finally disappear. 

In contrast to hard- wall billiard systems, the decay of 
the ray intensity in open optical system is an additional 
important characteristics. When the ray hits the dielec- 
tric boundary, the ray intensity decays according to Fres- 
nel's equation [3l|. The decay rate 7 is defined by 



-it 



(7) 



where I(t) is the remaining ray intensity after time t 
which is also a measure of the trajectory length. In order 
to obtain the decay rate 7 of attractor in RMDS, we con- 
sider an ensemble of initial points distributed uniformly 
over the open region of phase space. After the transient 
time, we calculate the remaining ray intensity accord- 
ing to Fresncl's equation whenever the discarded ray in 
RMDS emits from the system. The time t is scaled to 
be the length of ray trajectory inside two disks and 1/n 
of the length of ray trajectory outside disks. The decay 
rates of regular islands of Fig. [2] are found to be larger 
than 0.46 and smaller than 0.52 from numerical calcu- 
lations. In consistency with these values we obtain the 
decay rate ^hbb of the horizontal bouncing ball type 
periodic orbit as 0.4883 from Eq. ([7]). 

Up to now we have studied the ray dynamics of two 
coupled dielectric disks employing a ray model with a 
specially chosen selection rule. We have found the for- 
mation of attractors as an essential feature. The ques- 
tion arises to what extent the RMDS model accurately 
describes the reality. To this end we now study the wave 
dynamics of the system and, following the concept of 
ray-wave correspondence, will pay special attention to 
the possibility of modes that correspond to the attractor 
structures of the ray model description. 



III. RAY- WAVE CORRESPONDENCE IN 
COUPLED OPTICAL MICRODISKS 

Our next interest is thus 'How do the classical struc- 
tures emerging from RMDS in coupled optical microdisks 
relate to resonance modes?' In a single optical microdisk, 
the properties of the imaginary parts of the complex res- 
onances are well explained by semiclassical analysis |43| . 
Resonance modes can be classified into two groups ac- 
cording to their loss which is represented by the imagi- 
nary part. Modes of the first group are high-Q WGMs 
with dominant angular momenta above the critical p c for 
total internal reflection and, consequently, only tunnel- 
ing leakage to the exterior. The second group consists of 
low-Q modes with dominant angular momenta below p c 
implying refractive leakage. Figure [5] (a) and (b) shows 
resonances with transverse magnetic (TM) polarization 



(a) 




FIG. 5: (color online) Resonances in (a) single and (c) cou- 
pled optical microdisks when n = 2.0 and d = 1.0 on complex 
plane. In (c), the black circles, red rectangles, green dia- 
monds, and blue triangles are resonances of which symmetric 
classes are EE, OE, EO, and OO, respectively. The distribu- 
tions P[Im(fc)] of the imaginary parts of complex resonances 
(b) in a single optical microdisk and (d) with EE-parity in 
coupled optical microdisks. The red (dark) bars represent 
high-Q resonances. 



of a single disk with n = 2 in the complex plane that 
are obtained from the Hclmholtz equation [441 ]. and the 
corresponding distribution of their imaginary parts. The 
complex wave number k inside the microdisk is given as 
nkoRji, where kg is the wave number outside. In this 
case, the resonances near Im(fc) <~ are very high-Q res- 
onances and form the first group indicated by the long 
red bar in Fig.[5](b). The resonances near Im(fc) ~ —0.55 
are almost bouncing ball type resonances with very low-Q 
factor. The distribution of these low-Q resonances rep- 
resenting the second group as well as the lower limit of 
imaginary values can be explained by semiclassical anal- 
ysis. 

For the coupled optical microdisk system, we obtained 
resonances belonging to four symmetry classes (EE, OE, 
EO, and OO with even (E) and odd (O) symmetries 
with respect to the vertical and horizontal symmetry 
axes of the two-disk system) using the boundary ele- 
ment method [45| . They arc shown in Fig. [5] (c) for 
the parameters d — 1.0, R = Rl/Rr = 1.0, and 
n = ul/uq = tir/hq = 2.0. The distribution of imag- 
inary parts of the resonances with EE-parity is presented 
in Fig. [5] (d) [4(| . In contrast to the two peak structure 
found for the single disk, there are now three prominent 
peaks in the distribution of imaginary parts of resonances 
in coupled optical microdisks. 
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First, the very high-Q WGMs with imaginary parts 
larger than —0.02 are almost the same as those in a single 
optical microdisk because the intcrdisk distance which is 
inversely proportional to coupling strength is sufficiently 
large. For the high-Q WGMs, the two-disk system can be 
considered as consisting of two individual microdisks and 
a coupling term. The coupling between two microdisks 
can be explained by optical tunneling or frustrated total 
internal reflection [47], The tunneling process oc- 

curs at the interdisk spacing which is a small part of 
the circular boundary and does not significantly influ- 
ence the resonance position. Also, the high-Q resonances 
with imaginary parts between -0.2 and -0.02 do not dif- 
fer a lot from those in a single optical microdisk but the 
resonances split according to the symmetry classes and 
these splittings tend to increase as the imaginary parts 
become smaller. It is natural that the coupling strength, 
represented by the amount of splitting, is proportional 
to the resonance loss because larger emission from one 
of the microdisks induces a larger coupling to the other 
microdisk. 

Secondly, the peak in the distribution near Im(fc) ~ 
—0.55 corresponds to the bouncing ball type resonances 
in a single optical microdisk and is therefore directly re- 
lated to the single disk result in Fig. [5] (b). Although 
the results presented are for the TM case, there is gener- 
ally no lower limit of imaginary values because the mode 
coupling to the environment is not intrinsically limited. 
Rather, the co upling between internal and external (or 
outer or shape p9M52| ) resonances, i.e. modes localized 
inside and outside the cavity respectively, causes as broad 
distribution of imaginary parts. 

The striking contrast between the distributions in sin- 
gle and coupled optical microdisks is the peak near 
Im(fc) ~ —0.25 in Fig. [5] (d). Those resonances do not 
have a partner in the single optical microdisk. In fact, 
these resonances originate from the new classical struc- 
ture obtained within the RMDS ray model, namely the 
attractors discussed above. This is confirmed in Fig. [5] 
that shows the near field intensity patterns of two typi- 
cal regular modes corresponding to the island structure 
(RMCI) of the RMDS attractors in Fig. [2] The patterns 
of Fig. [6] (a) and (b) resemble the horizontal bouncing- 
ball-type periodic orbit in Fig. Q] (c) and the ray tra- 
jectories of the RMDS in Fig. [51 respectively. Figure [5] 
(c) and (d) represent the generalized Husimi functions 
[531 ] of the mode in Fig. [6] (b) and show incident and 
emerging intensities, respectively, at the inner boundary 
of upper semicircle of the left disk. The detail shown 
can be directly compared to the PSOS shown in Fig. [2] 
(a), and we find indeed a convincing correspondence to 
the island structure shown there. The Husimi functions 
follow it very closely, except for the missing of the upper 
part of the right island in panel (c), indicating that there 
is (almost) no incident intensity coming from inside the 
left disk, that is actually accompanied by a difference 
in intensity in the upper and lower lobes of this island 
in panel (d) for the emerging Husimi function. As the 
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FIG. 6: (color online) Near field intensity patterns of reso- 
nance modes localized on regular islands with wave numbers 
(a) k = 45.0025 - i0.2451 and (b) k = 50.9680 - i0.2582. 
Black-red(gray)-yellow(light gray)-white colors indicate high 
to low intensity, (c) Incident and (d) emerging generalized 
Husimi functions taken at the inner boundary of the left disk, 
the detail shown corresponds to the PSOS representation in 
Fig. [2] (a). The dashed lines \p\ — p c = 1/n — 0.5 mark the 
critical momemtum for the onset of total internal reflection. 
The arrows of insets are the ray trajectories corresponding to 
positions, (s,p) ~ (0.8, ±0.3), of upper and lower lobes of the 
right islands. 



RMDS-based PSOS does, unlike the Husimi functions, 
not distinguish between incoming and outgoing rays it is 
in fact the sum of two Husimi functions that has to be 
compared to the PSOS island structure in Fig. [5] (a), re- 
sulting in a nice agreement. The ray trajectories of insets 
in Fig. [S] (c) and (d) show the correspondence between 
the near field pattern of Fig. [6] (b) and the Husimi func- 
tions. The solid- and dashed-line arrows represent high 
and low intensities, respectively. The convincing agree- 
ment found in phase and real space, and especially the 
clear peak near Im(fc) ~ —0.25 in Fig. [5](d) suggest that 
indeed many resonances are closely associated with the 
RMDS island structure, which justifies in turn the use 
and the validity of the RMDS that we introduced and 
applied in the previous Section. 

In order to elucidate the correspondence between the 
decay rate 7 of the islands and the loss of resonances, we 
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FIG. 7: (color online) The distributions of the imaginary parts 
of resonances in coupled optical microdisks when (a) n = 1.8 
and (b) n = 1.5 with d = 1.0. Insets are near field inten- 
sity patterns of typical modes corresponding to the classi- 
cal structure of RMDS, of which k = 45.7196 - iO.2084 and 
k = 47.9012 - i0.5649, respectively. The blue line is the dis- 
tribution numerically obtained from the decay rates of the 
underlying classical strucuture of Fig. [3] (a) and Eq. ((9|. 



introduce the intensity of resonance decay given by 



J(t) = 



; -i(Re(fc)+iIm(fe))t 



_ e 2Im(fc)i 



(8) 



where we take the speed of light c = 1 in vacuum without 
loss of generality. From Eq. (|7J) and Eq. (jHJ, we obtain 
the relation 



Im(fc) = - 7 /2. 



(9) 



Since the island structure has decay rates 7 between 0.46 
and 0.52 as we obtained in previous section, RMCIs have 
imaginary values between —0.23 and —0.26, in full agree- 
ment with the numerical result Im(fc) ~ —0.25 read off 
from Fig. [5](d). We conclude that the new classical struc- 
tures found within the RMDS, such as the regular islands 
of Fig. [2 play an important role and can host RMCIs. 
As a consequence, the distribution of the resonance losses 
is changed with respect to the uncoupled system. More- 
over, the RMDS explains well the losses of RMCIs as well 
as the near field intensity patterns. 

Figure [7J (a) and (b) show the distributions of imag- 
inary parts of complex resonances when n = 1.8 and 
n = 1.5, respectively. The distributions have also a three 
peak structure, but the central peak is less pronounced 
and broader than in the case of n = 2.0. The broad distri- 
butions originates from the chaotic-to-regular transition 
of classical attractor island structure when the refractive 



index n is increased from 1.5 to 1.8, cf. Fig. |3l Below 
n ~ 1.732, the PSOS in the open region < p c ) is fully 
chaotic. When n = 1.8, only small regular islands in the 
vicinity of the hexagonal shaped periodic orbit of Fig. [T] 
(d) appear in the PSOS. The inset of Fig. [7] (a) shows 
a corresponding RMCI (Im(fc) = —0.2084, this value is 
consistent with the decay rate 7 ~ 0.4 of the hexagonal 
shaped periodic orbit). For n = 1.5, the near field inten- 
sity pattern shown in the inset of Fig. [7] (b) is typical for a 
chaotic mode with a rather low Im(fc) = —0.5649 and cor- 
respondingly high intensity outside the two microdisks. 
The distribution of the imaginary parts of resonances in 
a fully chaotic system is associated with the decay rate 
of the underlying classical structure of Fig. [3] (a) , which 
is represented by the blue line in Fig. [7J (b). Recently, it 
has been investigated that this distribution can also be 
related to the escape rate for the phase space points near 
the chaotic repeller including Fresnel's laws [EH, (5f| . 

We now consider the dependence of the resonances and 
the loss distribution on the system parameters. For the 
case of the high-Q WGMs, the real parts of the reso- 
nances split as the interdisk distance decreases if the two 
microdisks are sufficiently close Q and is inversely pro- 
portional to the radius of a microdisk or independent of 
the radius in coupled nonidentical optical microdisks be- 
cause the dominant WGM locates on only one microdisk 
[i~3l ]. Concerning the horizontal bouncing ball type peri- 
odic orbit of Fig. [T] Fig. [8] (a) and (b) show the real parts 
of the corresponding RMCIs as function of the interdisk 
distance and the ratio of radii R of the two microdisks, 
respectively. Unlike the real parts of high-Q resonances, 
the real parts of RMCIs increase as the interdisk distance 
or/and the radius of the left microdisk decreases as shown 
in Fig. [5] (a) and (b), respectively. The black and red full 
lines in Fig. [5] (a) represent RMCIs which localize on the 
horizontal periodic orbit. They undergo an avoided res- 
onance crossing near d = 0.6: For lower d, the resonance 
is of type B and follows the red line; after the avoided 
crossing, it follows the black line and changes its char- 
acter to A-type. As the near field intensity patterns of 
the A- and B-modes in Fig.[H]show, both modes have the 
mode index I = 15 on the periodic orbit. The resonances 
depending on the radius of the left microdisk in Fig. [8] 
(b) are also RMCIs with localized intensity patterns on 
the periodic orbit; their near field intensity patterns are 
shown in Fig. [5] 

Using the quantization rule on the periodic orbit, we 
obtain the real part Re(fc) of a resonance as a function 
of the radii, Rl and Rr, of the microdisks, the interdisk 
distance d, and assuming closed boundary conditions at 
both ends of the periodic orbit, as 

(10) 



2{R L + R B ) + d/n 

with a = 21 and 21 — 1 for Dirichlet and Neumann bound- 
ary conditions, respectively. For TM polarization, the 
Re(fc) satisfies the relation 



< Re(fc) < fee, 
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FIG. 8: (color online) The real parts of resonances as a func- 
tion of (a) the interdisk distance in coupled identical optical 
microdisks and (b) the ratio of radii in coupled nonidentical 
optical microdisks when d — 0.1. The four near field intensity 
patterns are for A-, B-, C-, and D-modes which are repre- 
sented by black circles in (a) and (b). 



where fcjy and kjy are the wave numbers for Neumann and 
Dirichlet boundary conditions, respectively. In Fig. [5J the 
green dashed and blue dotted lines represent k rj and k n , 
respectively, illustrating that the relation Eq. (fTTj) is well 
satisfied. 

We finally mention that coupling-originated classical 
structures such as the attractor islands discussed here 
may also affect the near field pattern of high-Q modes. 
Whereas it is generally accepted (l6l - [2(ij that the low- 
intensity tail structure of the near field patterns of high- 
Q modes in chaotic microcavities are related to the 
chaotic repeller structure which corresponds to the un- 
stable manifold structure near the critical line, this might 
change when new regular structures appear in the phase 
space due to the presence of a second disk to which the 
light can be coupled. Then, a hybridization [56j of high- 
Q WGM (single disk) and low-Q RMCI (coupled disks) 
is possible and was confirmed in numerical calculations. 



In this case, the near field pattern of the high-Q mode is 
not determined by the chaotic repeller structure near the 
critical lines, but rather by elements of the low-Q mode 
leakage which are typically positioned in the center of the 
leaky region instead. 



IV. SUMMARY 

We have proposed the ray model with deterministic 
selection rule (RMDS) for coupled optical microdisks as 
a theoretical model to deal with the inherent ray splitting 
dynamics in optical systems. We found that new classical 
structures in the form of regular islands (representing the 
attractor of the RMDS model) occurred, for a certain 
range of system parameters, in phase space. We have 
elucidated the stabilities and the decay rates of these 
structures which depend on the system parameters such 
as interdisk distance, ratio of radii, and refractive indices 
of microdisks. 

We have confirmed the physical significance of struc- 
tures emerging from the RMDS by identifying resonances 
associated with the classical structures and investigating 
the ray- wave correspondence. The near field intensity 
patterns of the resonances resemble the trajectories of 
regular islands and the real and imaginary parts can be 
explained by quantization rules and decay rates of ray 
trajectories, respectively. We expect that the ray model 
with a specific deterministic selection rule will be useful 
for the description of resonance modes in other coupled 
systems and generally in complex physical systems with 
splitting dynamics. 
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